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Abstract 

Precision tests of the CKM mechanism and searches for new physics in the 
flavour sector require dedicated QCD calculations of decay widths and spec- 
tra. Significant progress has been achieved in recent years in computing in- 
clusive B decay spectra into light energetic partons. 1 briefly review dif- 
ferent theoretical approaches to this problem focusing on the determination 
of | V^t, | from inclusive semileptonic decays and show that this determination 
is robust. The largest uncertainty is associated with the value of the b quark 
mass. Finally I present new numerical results in the DGE resummation-based 
approach, now including O((ioa 2 s ) corrections. The results are presented for 
all relevant experimental cuts, from which a preliminary average is derived 
\V uh \ = (4.30±0.16(exp)t[J;^(th)l^(m 6 ))-10- 3 , where the PDG value of the 
b quark mass, m^ s = 4.20 ± 0.07 GeV, is assumed. 



1 Introduction 



Low-energy precision measurements, in particular precision determination of 
the CKM parameters and the branching fractions of rare decays, provide many 
valuable tests of the Standard Model. The resulting constraints on new physics 
are highly complementary to the direct searches at hadron colliders, and are 
expected to continue being so throughout the LHC era [1]. 

The experimental effort over the past few years by the B factories and the 
Tevatron has established the fact that CKM is the main mechanism of flavour 
and CP violation in the quark sector. This is now a field of precision physics. 
Measuring deviations from the Standard Model and further strengthening the 
constraints on new physics will require a continuous experimental effort [2-5] 
alongside corresponding progress on the theory side. 

The obvious example of the progress that was made is the precise mea- 
surement of the small angle (5 of the Unitarity Triangle, and its comparison with 
the short side of the triangle, \V u b/V c b\- While the former is directly sensitive 
to potential CP-violation beyond the Standard Model and can be measured 
experimentally with high precision without any theoretical input, the latter, be- 
ing determined by tree-level Weak (semileptonic) decays, is insensitive to new 
physics, however it heavily relies on theoretical calculations in QCD. Currently 
the measured sin(2/3) is not entirely consistent with \V u b/V c b\, introducing some 
tension into global fits [6,7]. This comparison is a crucial element in the big 
picture. 

Experimentally, both \V u b\ and \V c b\ can be measured either using an ex- 
clusive hadronic final state, e.g. B — » D*W and B — > nlD or by considering 
the inclusive rate, summing over all hadronic final states subject to some kinc- 
matical constraints. The two approaches involve different experimental and 
theoretical tools and are therefore complementary. 

The lower rate of the b — > u transition (iV^/Fcbl 2 ~ 1/50) makes the \ V u b\ 
measurement more challenging. Further difficulties, common to all heavy- 
to-light decays, lie on the theory side. The exclusive determination of \V u b\ 
requires a theoretical calculation of the form factor using non-perturbative 
methods such as Lattice QCD [8,9] or QCD sum-rules [10], which both have 
systematic uncertainties that are hard to quantify. The Lattice calculations are 
expected to improve in the future. In contrast the inclusive determination relies 
primarily on the Heavy Quark Expansion (HQE) and QCD perturbation theory, 
where a systematic improvement can be achieved and uncertainties are easier 
to quantify. However, the case of inclusive b — > u decay is further complicated 
by the fact that the final state is characterized by jet-like kinematics. This, 
in conjunction with the experimental requirement to perform the measurement 
subject to stringent kinematic cuts (in order to suppress the charm background) 
implies that to extract \V u b\ one needs a precise calculation of the spectrum, 
not just the total width. This provides one of the biggest challenges in Heavy 



Flavour physics, a subject an intense theoretical effort over the past few years. 
This is the subject of the present talk. 

2 A brief look at inclusive semileptonic b — > c decay 

In order to appreciate the difficulty in determining \ V u b\, it is useful to compare 
it to the much favored case of \ V c b\- Let us therefore briefly review the situation 
in inclusive semileptonic b — > c decays. 

Owing to the high rate of these decays and their distinct characteristics, 
the B factories provide precise measurements of the branching fraction, as well 
as the first few moments, over the entire phase space. These truly-inclusive 
observables can be readily computed using the HQE, for example, 



T(B -» XJi>) = 




on — shell b— quark decay with IR cutoff 

where the first term stands for the partonic on-shell &-quark decay width T(b — > 
X C W), computed with an infrared cutoff fi, and other terms, suppressed by 
powers of the b— quark mass correspond to matrix elements of local operators, 
computed with an ultraviolet cutoff //,; as usual, the \i dependence cancels out 
order by order. Importantly, non-perturbative corrections first appear at order 
1 /m\, where there are two non-perturbative matrix elements, the kinetic energy 
H 2 . and the chromomagnetic energy hq. 

Note that the partonic decay width computed to next-to-leading order 
(order a s ), without a cutoff, already yields a viable approximation to the 
measured width. This approximation is systematically improved by including 
non-perturbative corrections as well as higher-order perturbative corrections. 
There is a continuous progress on both these fronts. Recent fits to the measured 
moments [11-13] are based on 0(/3oa 2 ) accuracy [14] with power corrections 
through 0(1/ ml), computed with leading-order coefficient functions. These 
fits (performed in the "kinetic" or IS mass schemes) provide a determination 
of \V c b\, mb and m c at the 1 — 2% level together with a determination of fi 2 at 
the 10% level. 

Very recently, complete 0(a 2 s ) corrections have become available in a fully 
differential form [15,16]; 0(1/ mf) corrections have been computed for the first 
time [17] and also the 0(a s ) correction to the coefficient function of [i 2 [18] 
was calculated. These advances will further push the accuracy, exploiting the 
potential of the B factory measurements for inclusive b — > c decays. 



3 The challenge: computing the charmless semileptonic decay spec- 
trum 

Experimental measurements [19-25] of inclusive b — > u decays involve kine- 
matic cuts in order to remove the charm background. Therefore, extracting 
\V u b\ from the data requires theoretical predictions for the fully (triple) differ- 
ential B — > X U W spectrum. 

The main difficulty in computing the spectra of heavy-to-light decays, 
B — > X u lv or B — > X s ~f, is in the fact that most events have a jet-like final 
state where the hadronic system X has a mass which is much smaller than 
its energy (approximately half of the energy released in the decay, mi,). The 
jet kinematics is most easily described in terms of light-cone momenta P ± = 
Ex TPx, where a typical event has P~ of order m& and P + not far above the 
QCD scale A. The decay process involves dynamics on scales that are far apart 
P~ 3> P + , complicating the perturbative description as well as the separation 
and parametrization of non-perturbative effects. 

It has long been recognized [26,27] that an attempt to compute the spec- 
trum in the small-P + limit by means of the HQE would run into serious diffi- 
culties: the dynamics is dominated by gluons with momenta of 0(P + ), turning 
the A/mft expansion into a A/P + one! The physical picture behind this break- 
down of the expansion is clear: the small lightconc component of the jet is 
influenced by soft gluon radiation as well as small fluctuations in the momen- 
tum carried by the decaying heavy quark. To recover a useful heavy-quark 
expansion, the dominant effects, those controlled by the scale P + , must be re- 
summed to all orders. This sum gives rise to the well-known "shape function", 
which can be interpreted as the momentum distribution function of the b quark 
in the B meson. Similar functions appear at higher orders in the heavy-quark 
expansion. 

Recall that in the b — > c case one could make use of the HQE: the rele- 
vant observables were well approximated by perturbation theory and the non- 
perturbative corrections were restricted to a few local matrix elements. In 
contrast, when considering the b — > u case one is required to compute the 
spectrum at small P + , which is proportional, already at leading power, to a 
non-perturbative object, the "shape function". This function is defined by the 
non-local matrix element: 

S(k+; n) = J™ ^e"^" (B\ h(y)[y, 0} 1+ h(0) \B) , (2) 

where k + is a lightcone momentum component (h are heavy-quark effective 
theory fields, and [y, 0] represents a gauge link). This function describes the 
distribution of momentum carried by the b quark. Thus, we observe that 
instead of having a few unknown non-perturbative matrix elements which enter 



as power-suppressed corrections, one faces here an unknown function already 
at the leading order in A/mj,! 

Described in these terms the problem of computing the spectrum and 
extracting \V u b\ from data may appear hopeless, or at least require a full- 
fledged non-perturbative approach. In fact, as we shall see below, the actual 
situation is significantly better. The partial branching fractions corresponding 
to experimentally relevant cuts, which vary between 20 to 60 percent of the 
total, can be still estimated reliably with very little non-perturbative input. 
Moreover, at present, the largest uncertainty in extracting \V u b\ is associated 
with the parametric dependence on the b-quark mass; other uncertainties (e.g. 
power corrections, Weak Annihilation) can be reduced by further exploiting 
the data and thus the prospects for an even more precise \V u b\ from inclusive 
decays are high. 

In the following I will briefly describe different theoretical approaches that 
have been developed in the past few years to compute the fully-differential 
spectrum and thus extract \V u b\ from the B factories data. I will not enter 
into any technical details, just try to give the flavor of the physics involved 
and the principal differences between the approaches. I will also not cover 
all the interesting theoretical developments in this area, notably the method to 
express the b — > u branching fractions directly in terms of the measured photon- 
energy spectrum in B — > X s -f, which had some resurrection recently [28-30], 
incorporating subleading effects in A/mj,. 



4 HQE based structure function parametrization approach 

The central idea of the HQE-based structure-function parametrization ap- 
proach, which has been recently put forward and implemented by Gambino 
et. al. [31], is to first use the HQE to compute carefully-selected observablcs 
- the first few moments of the structure functions — where this expansion 
is expected to be most reliable, and then use these observables to constrain 
the parametrization of the spectrum. In this way one bypasses the need to 
deal with the difficult kinematic region P~ ^> P + where neither the HQE nor 
perturbation theory converge well. 

To explain briefly how the calculation in this approach is set up, let us 
recall that the triple-differential rate can be written in terms on three hadronic 
structure functions Wi(qo, q 2 ): 



dT G 2 F \V U ; 12 



dq 2 dq dE l 8tt 3 



2E 2 - 2q Q E l + Y 



W 2 + q 2 (2Ei - q )W 3 



(3) 

where qo and q 2 are the total leptonic energy and squared invariant mass, 
respectively. Ref. [31] computes the shape of the physical structure functions 
Wi(qo,q 2 ) as a convolution at fixed q 2 between non-perturbative distribution 



functions Fi(k + , q 2 ; /j,) and the perturbative (presently the Born- level) structure 
functions W? eit {q Q ,q 2 ): 



W t (q ,q 2 ) = J dk+F^q^-nWr' (<*> - ^ (l - ,q 2 ;^ , (4) 

where the functions Fi(k+, q 2 ; /i) are parametrized and constrained by the first 
few go-moments of Wi(qo, q 2 ). 

The moments used to set these constraints are computed using the HQE, 
where the perturbative part currently includes corrections up to O(a 2 /3o) [32] 
and power corrections are included through 0(1/ ml). The separation between 
the perturbative component and the power-correction terms is based on a 
hard momentum cutoff (/j, = 1 GeV) in the "kinetic scheme", which has the 
advantage that the input parameters (in particular and /ij) can be taken 
directly from fits to the b — > c moments. 

Similarly to the b — > c analysis both the power expansion and the pertur- 
bative expansion^ can be improved once higher-order corrections are known. 
The use of a hard cutoff on the gluon energy (the "kinetic scheme" ) eliminates 
the sensitivity to multiple soft emission rendering the expansion better conver- 
gent. Nevertheless a single-logarithmic collinear divergence persists, and can 
in principle be resummed. 

In this approach the parameters in Fi, i — 1,2,3 are fixed a new at each 
given value of q 2 , based on the moment constraints. Thus, the way these pa- 
rameters vary with q 2 is indirectly determined by the HQE. This issue becomes 
crucial at large q 2 , where the HQE breaks down: the final-state hadronic sys- 
tem is then soft. While the contribution from this phase-space region is small 
(it is power suppressed) the spectrum there is clearly not well under control. 
Ref. [31] provides an interesting analysis of the breakdown of the HQE in this 
region and relates it to the presence of 0(1/ ml) Weak Annihilation contribu- 
tions, centered at q 2 ~ m 2 . The contributions from the large-g 2 region are 
parametrized making a conservative estimate of their size. Even then the im- 
pact of the Weak Annihilation contributions on the average value of \ V u b\ is just 
~ 3% [33], which is less than the parametric uncertainty due to m-t dependence. 
Having said that, Ref. [31] has clearly demonstrated that further experimental 
input on the q 2 distribution and moments, measured separately for charged and 
neutral B mesons, would be important for reducing the uncertainty on |V^b|. 

To summarize, the approach of Ref. [31] is cautious: it uses the well- 
understood (and well tested!) theoretical framework of the HQE for carefully 
selected moments, and assumes very little beyond that. It relies however, on 

x At present fixed-order O(a 2 {3o) expressions are used. It is fair to say that 
the generalization of the hard cutoff approach beyond the level of a single gluon 
(possibly dressed) is difficult to implement. 



extensive parameterization, dealing with three non-perturbative functions Fi 
of two kinematic variables, whose properties are unknown. The authors of 
Ref. [31] therefore took special care to consider a large class of functions and 
further devised means to assess whether this class is large enough. In this way 
they managed to provide a reliable prediction for the triple differential spectrum 
over the entire phase space without dealing directly with the difficult kinematic 
region where the hadronic system is jet-like. In the following I present other 
theoretical approaches that instead consider directly this region. 

5 Shape function approach 

The shape-function approach by Neubert and collaborators [34, 35] deals di- 
rectly with the important kinematic region where P + <C P~ ■ This is done by 
establishing a modified expansion in inverse powers of the mass, where at each 
order the dynamical effects that are associated with soft gluons, fc + ~ 0(P + ) ~ 
0(A), are summed into non-perturbative shape functions. At leading power 
there is one such function, the momentum distribution function defined in ^ 
above; beyond this order there are several different functions with additional 
fields insertions. To extend the calculation beyond this particular region, it is 
constructed to match the standard HQE when integrated over a significant part 
of the phase space. In this way two systematic expansions in inverse powers of 
the mass are used together. 

The modified expansion in shape functions is developed, following the 
Soft Collinear Effective Theory (SCET) methodology [36-38], for the partic- 
ular kinematic region where the final-state is jet-like P~ ~ m&, P + — A, 
and thus mx = \/ TO &A, the region into which a large fraction of the events fall. 
The large (parametric) hierarchy between these scales implies loss of quantum- 
mechanical coherence between the respective excitations, leading to factoriza- 
tion into three different subprocesses [39] (see also [34,36,40]). The result, at 
leading power in A/mb, can be expressed as a convolution integral [35]: 



where y = (P~ — P + )/(Mb — P + )- Here H stand for hard, depending on 
momenta of 0(m&), J for jet, depending on momenta of C(\/ m ftA) and S for 
soft, depending on momenta of 0(A) and is therefore considered as a non- 
perturbative object, to be parametrized. Similar factorization formulae apply 
at subleading powers in A /mj, leading to the following expression for the dif- 
ferential width: 



dT 



p+ 



dP-dP+dEi 



dk + ym b J (ym b (P+ - k + ),{i) S(k+, fi) (5) 



(6) 



The matching into the standard HQE translates into constraints on the mo- 
ments of the shape functions. 

The authors of Ref. [35] have defined the separation between the leading 
term and the power corrections using a factorization procedure that is based 
on dimensional regularization. They parametrize the shape functions directlj0 
at the intermediate (jet) scale [i ~ \frn^K (in practice /i = 1.5 GeV is used) 
and thus avoid ever dealing with softer momentum scales. 

Both the hard and the jet functions are computed in perturbation the- 
ory. Owing to the presence of well-separated scales, the jet energy 0(mt>) and 
its mass 0(^fm^K) the perturbative expansion contains large Sudakov loga- 
rithms. Ref. [35] resum these logarithms to all orders with high logarithmic 
accuracy (next -to-next-to-leading log, NNLL). The hard coefficient function 
is currently computed at 0{a s ). 

Variation of the matching scale is used to estimate missing higher-order 
corrections. This translates into about 3 — 4% uncertainty on the average value 

of \V ub \. 

The functional forms of the shape functions are unknown, and they are 
therefore parametrized. The first two moments of the leading shape function are 
reasonably well constrained: they are fixed by toj, and fi^ respectively; higher 
moments are not well constrained. Subleading shape functions are difficult to 
constrain. 

To summarize, the method of Ref. [35] makes extensive use of the avail- 
able theoretical tools, employing consistently two expansions that are valid in 
two different kinematic regimes: the expansion in shape functions, valid for 
the typical final-state momentum configuration, and the standard HQE, valid 
for the fully integrated width. Further to that, Sudakov resummation for the 
jet-scale logarithms is employed at NNLL accuracy. By using a relatively high 
factorization scale the perturbative calculation remains insensitive to the in- 
frared, and converges well. All ingredients, the power expansions as well as the 
perturbative ones, can in principle be improved systematically by including 
higher-order terms. 

The motivation to go to subleading powers in a completely general way, 
however, has a price: one needs to parametrize several different subleading 
shape functions on which there is no theoretical control nor experimental 
input. Despite this, the authors of Ref. [35] have demonstrated that the 
experimentally-relevant partial branching fractions remain under good con- 
trol: their sensitivity to the unknown higher moments of the leading shape 
function as well as to the unknown functional form of the subleading shape 
functions is small: the estimated effect of these unknowns on the average value 
of \V u b\ is less than 1%! The largest uncertainty in the determination of \V u b\ 

2 Note that in this way the (known) evolution properties of the shape- 
function are not being used. 



is the parametric one, owing to the strong dependence on mi,. 

A common feature of the two approaches described so far is the extensive 
use of parametrization of non-perturbative functions. The fact that these func- 
tions all have a clear field-theoretic definition does not presently help in their 
parametrization. To improve on that one needs to avoid introducing an explicit 
cutoff nit, /i ^> A. This is indeed possible. It is well-known that inclusive 
observables, such as moments of decay spectra, are infrared-safe observables. 
In other words, in the absence a cutoff soft-gluons divergences cancels out 
in the sum of real and virtual diagrams, making the moments finite at any 
order in perturbation theory. In the following I shall describe a resummation- 
based approach, where a cutoff is not used and consequently one relies less on 
parametrization. 

6 Resummation based approach 

The approach of Refs. [40-43] uses resummed perturbation theory in moment- 
space to provide a perturbative calculation of the on-shell decay spectrum in the 
entire phase space without introducing any external momentum cutoff; non- 
perturbative effects are taken into account as power corrections in moment 
space. Resummation is applied to both the 'jet' function and the 'soft' (quark 
distribution) function, dealing directly with the double hierarchy of scales char- 
acterizing the decay process. Consequently, the shape of the spectrum in the 
kinematics region where the final state is jet-like is largely determined by a 
calculation, and less by parametrization. 

The resummation method employed, DGeJH, combines Sudakov and renor- 
malon resummation. Sudakov logarithms are resummed with high logarithmic 
accuracy (NNLL) for botlfl the 'jet' and the 'soft' functions. 

Renormalon resummation is an essential element in implementing a con- 
sistent separation between perturbative and non-perturbative corrections at 
the power level. Refs. [41-43] have adopted the Principal Value procedure to 
regularized the Sudakov exponent and thus define the non-perturbative pa- 
rameters. This is in full analogy [49] with the way a momentum cutoff is 

3 Dressed Gluon Exponentiation (DGE) is a general resummation formalism 
for inclusive distributions near a kinematic threshold [44]. It goes beyond 
the standard Sudakov resummation framework by incorporating renormalon 
resummation in the calculation of the exponent. This has proven effective [47, 
48] in extending the range of applicability of perturbation theory nearer to 
threshold and in identifying the relevant non-perturbative corrections in a range 
of applications [44-48] . 

4 The 'jet' logarithms are similar to those resummed in the approach of 
Ref. [35]; there however 'soft' logarithms are not resummed. 



conventionally used. Most importantly, this definition applies to the would-be 
l/m& ambiguity of the 'soft' Sudakov factor, which cancels exactly [40] against 
the pole-mass renormalon when considering the spectrum in physical hadronic 
variables. The same regularization used in the Sudakov exponent must be ap- 
plied in the computation of the b-quark pole masfl This vital mechanism 
is absent in a fixed-logarithmic-accuracy procedure (as employed for example 
in [50,51]) leading to an uncontrolled shift of the entire spectrum with P + /P~ . 

Aiming to provide a good description of the spectrum in the kinematic 
region where there is a large hierarchy between the lightcone momentum com- 
ponents, P + <C P~ , it proves useful to consider the moments with respect to 
the lightcone-component ratio^|: 

dT N (p-, E t ) _ / P +\ N - 1 dT(p+ 1 p" Ei) 



dp dE\ Jo ^ \ P ) dp + dp dEi ' ^ ^ 

where the partonic lightcone momentum components p^ are related to the 
hadronic ones by: p ± = P ± — A, where A = Mb — nib is the energy of the 
light-degrees-of-freedom in the meson. 

Note that in ([7]) large moment index corresponds to the limit of interest, 
jet kinematics: the main contribution to the integral for N — ► oo comes from 
the region where p + jp~ — > 0. For large N one identifies three characteristic 
scales, hard 0(p~), jet 0(p~ /vN) and soft (D(p~ /N). In this limit, and 
up to \JN corrections, the moments factorize [34,36,39,40] to all orders as 
follow^: 

"' T A ''" &) - H(p-, E t ) J( P -/VN, (jl) S b (p-/N, n) +0(l/N), (8) 



dp dE[ 



Sud(p-,N) 



where the factorization-scale (fi) dependence cancels exactly in the product in 
the Sudakov factor Sud(p~, N). 

To use perturbation theory one must consider p~ large enough, such that 
even the 'soft' scale p~ /N is sufficiently large compared to the QCD scale 
A. This hierarchy of scales is illustrated in figure [TJ Here one should note a 



5 In Eq. (fT0|) below the cancellation of the renormalon ambiguity involves 
Sudakov factor of Eq. ([8]), on the one hand, and A on the other. 

6 Note that we consider here the moments of the fully differential width [42]: 
the moments remain differential with respect to the large lightcone component 
p~ as well as the lepton energy Ei . This is essential for performing soft gluon 
resummation. This issue has also been discussed in Ref. [51]. 

7 Note that this factorization formula maps directly onto eq. ([5]) above, where 
the convolution integral turns into a product in moment space. 



subtle but important distinction from the shape-function approach discussed 
above, where it was a priori assumed that the "soft" scale (here p^ /N) is 0(A), 
prohibiting any perturbative treatment of the corresponding dynamical subpro- 
cess. Here instead we wish to compute Sb(p~ /N, /i) — the quark distribution 
inside an on- shell heavy quark [52, 53] — in perturbation theory, as a basis 
for the description of the physical distribution, the quark distribution in the B 
meson Sb(p~ /N 1 fi). Because Sb(p~ /N, fi) is infrared safe Sb(p~ /N, fi) only 
differs from Sb(p~ /N 1 fij by power corrections, powers of NA/p~ . Eventually, 
at N 3> p~/A all these powers become relevant, recovering the "shape func- 
tion" scenario. Refs. [41-43] therefore parametrize these power corrections. It 
should be noted that experimentally-relevant branching fraction are not so sen- 
sitive to the high moments, and therefore the effect of these power corrections 
is small. The resulting effect on \V u b\ is at the sub-percent level. 
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Figure 1: The hierarchy of scales underlying factorization as conceived in 
a cutoff-based approach (left) vs. the moment-space resummation-based ap- 
proach (right). 



Factorization facilitates the rcsummation of Sudakov logarithms, the cor- 
rections that dominate the dynamics at large N [34,36,39,40,51,54,55]. There 
is, however, another class of large corrections which is always important at 
high orders: these arc running coupling corrections, or renormalons. In the ap- 
proach of Refs. [41-43] the Sudakov exponent is computed as a Borel sum, facil- 
itating simultaneous resummation of Sudakov logarithms and running-coupling 



corrections. The Sudakov factor takes the form: 



, , , ,C F f du / A 

Sud(p-, AT) = exp ^ — - / — — 

Po Jo u \P 
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where G(u, N) = T(-u) 



T(iV) 

r(N - u) r(i - u) 



B 5 (u)G(2u, AT) - Bj(u)G(u, N) 
1 



(9) 



Here B$(u) and Bj(u) are the Borel representations of the Sudakov anomalous 
dimensions of the quark distribution and the jet function, respectively. These 
functions are known [53] to NNLO, 0(u 2 ), facilitating Sudakov resummation 
with next -to-next-to-leading logarithmic accuracy [41]. 

Beyond that, the analytic structure of the integrand is indicative of power 
corrections. The Sudakov exponent has renormalon singularities at integer and 
half integer values of u, except where Bs,j{u) vanish. The corresponding am- 
biguities, whose magnitude is determined by the residues of the poles in ((5J), 
are enhanced at large N by powers of N. They indicate the presence of non- 
perturbative power corrections with a similar N dependence. These power 
corrections exponentiate together with the logarithms. By evaluating the Borel 
integral, rather than expanding it, Refs. [41-43] make use of this additional 
information, defining the perturbative part of the exponent via the Principal 
Value (PV) prescription, and then parametrizing the dominant power correc- 
tions. 

So far only non-perturbative corrections that are leading in the large N 
limit, (NA/p~) k for any k, have been taken into account in this approach. 
These power corrections are the non-perturbative content of the leading shape 
function in the approach of Sec. [5] 0(1/N) effects corresponding to sub-leading 
shape functions in {B} are only accounted for in Refs. [41-43] at the perturbative 
level. In principle, however, subleading non-perturbative effects can also be 
parametrized as power corrections in moment space. This would be worthwhile 
doing at the point where the constraints on the leading power terms would be 
sufficiently tight. 

Moment space proves convenient for resummation and parametrization of 
power corrections, but at the end of the day one needs the spectrum in momen- 
tum space. The fully differential spectrum in hadronic variables is obtained by 
an inverse Mellin transform (c/. (J7J): 
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( 10 ) 

where the integration contour runs parallel to the imaginary axis, to the right 
of the singularities of the integrand. 



To summarize, the moment-space resummation approach of Refs. [40-43] 
allows to compute the fully-differential spectrum in the entire phase space as 
an infrared-safe quantity, without introducing any explicit cutoff scale. Non- 
perturbative effects are treated as power corrections, where the parameters 
are defined using the Principal Value prescription. This approach thus max- 
imizes the predictive power of perturbation theory, and minimizes the role of 
parametrization. In the next section we shall have a quick look at the resulting 
phenomenology. In particular, we will present here for the first time numeri- 
cal results that are based on matching the resummation formula to Q((5qol 2 s ), 
incorporating the results of Ref . [32] . 

7 \V u b\ by DGE including O(/3 a 2 ) corrections 

So far the calculation of the partial branching fractions from which \V u b\ was 
extracted, has been based on a NLO result: although the jet and the soft 
functions were resummed with NNLL accuracy, the hard coefficient function 
H(p~, Ei) in ©, corresponding to constants and 1/N suppressed terms at 
large N, was only known to NLO, 0(ct s ) [56]. In a recent paper [32] we have 
computed analytically the running-coupling corrections, which are the domi- 
nant corrections at the NNLO. Both real and virtual O((3oa 2 ) corrections are 
now available. 

Very recently I have completed the task of matching the resummed triple- 
differential rate to the new O((3oa 2 ) corrections, and implemented it into the 
C++ DGE program. The details of the matching procedure will be published 
separately. The new version of the program (Version 2.0) is available at [57]. 
Preliminary results based on this new version will be presented below. 

Prior to describing the new results a comment is due concerning the way 
the partial branching fractions are computed, which has changed between the 
old and new implementations [57]. In the old version the triple differential 
width, normalized as 1 /T^dT / dP + dP~ dEi {Tq is the Born-level width) was 
integrated over the relevant phase-space, and then divided by a normalization 
factor corresponding to a similar integral over the entire phase space. In the 
new version, I apply the same procedure but this time evaluating at each point 
in phase space the (perturbatively) normalized rate l/T to teddT / dP + dP~ dEi 
instead of l/T^dT / dP + dP~ dEi. This implies that the expression for ro/r to tai 
has been expanded, and multiplied into the hard matching coefficient. Finally 
the new hard matching coefficient is truncated at the required order, a s (NLO) 
or Poot 2 (NNLO). When working at NLO this amounts to an 0(a 2 ) difference 
with respect to the previous calculation, which is not small numerically (it is 
comparable to the Poet 2 term in the matching coefficient, and has the oppo- 
site sign). The new formulation is theoretically favored as it leads to smaller 



renormalization-scale dependency. 

Using the 2007 PDG value for the short-distance b-quark mass, 



mf 3 = 4.20 ± 0.07 GeV (11) 
I have computed the normalized partial widths 

_ T(B -> XJ9; cut) 

1 total (.-D — * AutfJ 

for the specific cuts used by HFAG [33] to extract \V u b\ based on the measure- 
ments in [19-25]. The results are summarized in table [TJ Note that the lepton 
energy cut is sometimes applied in the T(45) frame rather than the B rest 
frame, involving a boost of j3 = 0.064; this is indicated in the table and taken 
into account in the calculation. 

As shown in table [1] a significant contribution to the uncertainty is due 
to the parametric dependence on nib for which we have taken the conservative 
range of (fTTj) . Other uncertainties we take into account (all summed up in 
quadrature in table [J) are: 

• Parametric uncertainty in the input value of a s , where we take a^ ls (Mz) = 
0.1176 ± 0.020. 

• Power corrections associated with the quark distribution function, esti- 
mated by varying the u = 3/2 renormalon residue as well as the power 
terms based on the parametrization presented in Sec. 4.3 in [43]. We 
take (C3/2, / PV ) = (1,0) as default and determine the uncertainty by 
considering the case (C 3 / 2 ,/ PV ) = (6.2,0.3). 

• Weak Annihilation effect. We assume that Weak Annihilation effects can 
increase the width by up to 2%. This error is taken as unidirectional. 

• The residual dependence on /i in the matching coefficient is used to esti- 
mate higher-order perturbative corrections. We vary it from fj, = nib/2 
to fi — 2rnb, where the central value is taken at fi = mt,. 

Note that the NNLO result for each of the cuts is consistent within errors 
with the NLO one. For NLO we only quote the central values; the errors are 
similar to those at NNLO. In particular, considering here the normalized i? C ut 
computed by integrating l/Y tota \dT / dP + dP~ dEi the renormalization-scale de- 
pendence is low (1-2% on \V u b\) already at NLO and there is no significant 
improvement going to NNLO. 



8 Renormalization scale dependence appears in our formulation only through 
the hard matching coefficients, as running coupling corrections are resummed 
in the jet and soft functions. 



Table 1 : Computed values of R cu t for different experimentally relevant cuts. 
The collumn on the left described the cuts. For R^t° we present the central 
value only, while for i?^ LO the errors are broken into ( an asymmetric ) theory 
error (which includes parametric uncertainty in the input value of ot s ; renor- 
malization scale uncertainty; power corrections associated with the quark dis- 
tribution function; and Weak Annihilation effect) and parametric uncertainty 
in the input value of mb according to ill]) . 



cut 


Ref. 


pNLO 
"cut 


pNNLO 
"cut 


Ei > 2.1 GeV 
£, T(4S) > 1 GcV; m x < 1.7 GeV; q 2 > 8 GeV 
£ ; T(4S) > 1.9 GeV 
£ ; T(4S) > 2.0 GeV 
Ei > 2.0 GeV; ^ iax < 3.5 GeV 2 
£, T(4S) > 1.0 GcV; m x < 1.7 GeV 
£, T(4S) > 1.0 GcV; m x < 1.55 GeV 


[19] 
[20] 
[21] 
[22] 
[23] 
[24] 
[25] 


0.234 
0.372 
0.389 
0.311 
0.239 
0.658 
0.559 


U.ZZD _ .013 -0.024l'' i 6J 
n Qfifi +0.021 +0.016/ \ 

n n 7/1 +0.021 +0.023/ \ 
0.374 _ 013 _ .022( TO b) 

n qm +0.020 +0.024/ s 

0-301 _ 015 _ .025( TO b) 

232 + 020 +0 023 (m h ) 

U.zoz _ .013 -O.Q22\" l b) 

R9S+ 026 +°- 055 ('rn^ 
U.OZS _ .017 -0.059 l TO 6 J 

^2 +0 - 028 + 069 (m,\ 
U.OOZ _ 023 -0.071 \' n b) 



Next, to extract |V^b| the experimental partial branching fractions [19-25] 
can be directly compared to the theoretical calculation: 

AB(B -> XJD; cut) = r B x r total (B -» XJv) x R cut . (13) 

Using the calculation of Sec. 2 in [42] with PDG value of mj (jTTJ) we get the 
following value for the total width: 

777^72 r to tai(S -> X U W) = 67.3 ± 5.4PS" 1 . (14) 

I Kb I 

Using the updated world average value of the B-meson life time, tb = 1.573 ps, 
together with (fT4]l and the i? cut values of table Q] we obtain for |V^f,| the values 
quoted in table [2l Note that the uncertainty in the total width is dominated by 



Table 2: Extracted values of \V u b\ based on the measured partial branching 
fractions (see quoted references), using Eq. U3\) with the R cu t values of tabled 
and the total width of Eq. |_?^[ ). The errors quoted for NNLO are experimental 
(statistic and systematic raised in quadrature); theoretical, through R cu t; and 
parametric dependence on nib in both the total width and through R cu t- 



cut 


Ref. 


\v ub \ 

NLO 


\v ub \ 

NNLO 


Ei > 2.1 GeV 

£, T(4S) > 1 GeV; m x < 1.7 GeV; q 2 > 8 GeV 

£ ; T(4S) > 1.9 GeV 
E r(4S) >2 Q GeV 

Ei > 2.0 GeV; S£ iax < 3.5 GeV 2 
E? iiS) > 1.0 GeV; m x < 1.7 GeV 
£, T(4S) > 1.0 GeV; m x < 1.55 GeV 


[19] 
[20] 
[21] 
[22] 
[23] 
[24] 
[25] 


3.64 
4.33 
4.54 
4.17 
4.18 
4.21 
4.47 


3.73±0.44(exp)tg;iJ(th)lg;||(m 6 ) 
4.36±0.46(exp)^(th)i°;^(m 6 ) 

4.63±0.43(exp)tH2(th)^.3i(^) 
4.24±0.29(exp)l° : ^(th)l°;i(m 6 ) 
4.24±0.29(exp)t^ 2 7 (th)t«; 4 g(m b ) 
4.31±0.29(exp)t^(th)t[; : 4 «(m ;) ) 
4.58 ± 0.22(cxp)t|]^(th) 4 :0;l 4 (m b ) 



nib and it is therefore fully correlated with the parametric uncertainty associ- 
ated with nib in i? cu t . This is taken into account in the parametric uncertainty 
quoted in Tabled Because the effect of changing mj, on R cut acts in the same 
direction as in the total width, their product, which enters the determination of 
| Kb | in (fTB")) . is highly sensitive to m b . This is clearly reflected in the parametric 
uncertainty quoted in the table. It is also illustrated in figured 

Examining the values of \V u b\ corresponding to different cuts one observes 
very good agreement. Even ignoring the theoretical and parametric uncertain- 
ties (which are correlated), they all agree very well. There is one case where 
the agreement is not as striking: this is the CLEO result of Ref. [19] where 
the central value falls below all other determinations; note however the large 
experimental error quoted. 



Averaging the result^ in table [5] we obtain 

\V ub \ = (4.30 ± 0.16(fflq))lg;g(th)tg:»(m 6 )) • lO" 3 . (15) 

This result can be compared with other theoretical methods used to compute 
the partial widths (we only refer here to the two other methods that have 
been discussed above; HFAG [33] presents additional results). The HQE-based 
parametrization of Ref. [31] yields an average [33] 

\V ub \ = (3.94 ± 0.15(exp)t°;20(th)) ■ 1CT 3 , (16) 

while the shape-function approach of Refs. [34, 35] yields [58] 

\V ub \ = (4.31 ± 0.17(exp) ± 0.35(th)) • 1(T 3 , (17) 

where in all cases we quoted the numbers corresponding to the central value of 
Eq. (jlip — because of the strong m b dependence this requirement is essential for 
any valuable comparison. For the average value one finds very good agreement 
between dl~5l) and (|17[) and compatibility with (fl7J)) . We note that for the mx~ 
based cuts there is better agreement between the methods of Ref. [31] and 
Ref. [58], which both yield somewhat lower central values for |V^j,| (\V ub \ — 
4.0 • 10 -3 [33, 58] for the above m b ) as compared to the DGE approach. 

8 Conclusions 

I have given an overview of the main theoretical approaches used to compute 
the triple differential spectra in order to extract of \V ub \ from data. I have 
mainly emphasized the conceptual differences and the relations between the 
approaches, but I also reported briefly on their status, their formal accuracy 
and their particular sources of uncertainty. 

It is evident that despite making different approximations, the various 
determinations are consistent with each other. Add to that the remarkable 
consistency between different measurements that use different kinematics cuts 
- which provides a valuable confirmation for the theoretical description of 
the spectrum — the conclusion is clear: the inclusive determination of \V ub \ is 
robust. This puts us on firm grounds coming to examine the consistency of the 
CKM mechanism. 

Finally, the single most important source of uncertainty in the inclu- 
sive determination of \V ub \ is the b-quark mass. The dependence on the mass 

9 Note that in this average we have neglected several correlations. A more 
precise updated average is being prepared by the HFAG. 



is extremely high owing to the fact that both the total width and the cut- 
dependence increase with increasing mj. The effect this has on \V u b\ is shown 
in figure [5J Clearly, improving our knowledge of the b-quark mass would di- 
rectly translate into more precise \V u b\. 

Before concluding I find it appropriate to add a few words about the field. 
Beyond their obvious significance to phenomenology, inclusive B decays are a 
remarkable source of interesting theoretical problems in QCD. We have only 
scratched the surface of this exciting field in this talk. 

Inclusive decay are also very challenging experimentally, and although 
the experimental issues have not been mentioned here, it is obvious that there 
would not have been much point in giving this talk if not for the remarkable 
achievements of the B factories in this area. The on-going discussion between 
theory and experiment has also been extremely fruitful, and I would like to 
thank all those who have contributed to that. 
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Figure 2: The extracted value of \V u h\, averaged over different measurements 
with different kinematic cuts, as function of the b quark mass, m^ ls . The quark 
mass serves as an input to the calculation of the partial width; it affects the 
partial width thought the total width, ~ m^, and through R cu t, both acting in 
the same direction. The calculation of R cut is done by DGE including O(Poa^) 
NNLO corrections. The vertical error bars are based on the remaining sources 
of uncertainty, theoretical and experimental, added in quadrature. The wide 
red horizontal bar at the bottom shows the 1-sigma range of the PDG world 
average value setting the range of mt values we consider. Just above 

it we present two specific determinations of the mass: the one extending to 
the right (green) based on a HQE~based fit to inclusive moments of b — ► c 
decays [11, 12], converted to MS [58], and the other (blue) based on a recent 
precise determination [59] using the total cross section in e + e~ — > hadrons 
near the bottom production threshold. 



